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My lecture is devoted to the analytical results available for a large class of ax- 
isymmetric stationary flows in the vicinity of compact astrophysical objects. First, 
the most general case is formulated corresponding to the axisymmetric stationary 
MHD flow in the Kerr metric. Then, I discuss the hydrodynamical version of the 
Grad-Shafranov equation. Although not so well-known as the full MHD one, it al- 
lows us to clarify the nontrivial structure of the Grad-Shafranov approach as well 
as to discuss the simplest version of the 3 + 1-split language - the most convenient 
O |' one for the description of ideal flows in the vicinity of rotating black holes. Finally, 

I consider several examples that demonstrate how this approach can be used to 
obtain the quantitative description of the real transonic flows in the vicinity of 
rotating and moving black holes. 



Many astrophysical sources are axisymmetric and stationary to a good 
accuracy. These include both accreting neutron stars and black holes, ax- 



isymmetric stellar (solar) winds, jets from young stellar objects, and ejec- 
tion of particles from magnetospheres of rotating neutron stars. It can not 
be ruled out that such magnetohydrodynamic flows also play an important 
role in other galactic sources, e.g., microquasars. The latter are regarded 
as candidates for black holes not to say about active galactic nuclei where 
the electrodynamical processes in the vicinity of the rotating supermassive 
black holes are considered as the most reasonable model of their central en- 
gine 1,2 . So, it is not surprising that ideal magnetohydrodynamics, which 
allows sufficiently simple formalization of the problem, is actively applied 
when describing these flows. 

The point is that due to axial symmetry and stationarity (as well as the 
ideal freezing-in condition), in the most general case it is possible to intro- 
duce five integral of motions which are constant at axisymmetric magnetic 
surfaces. This remarkable fact allows us to separate the problem of finding 
the poloidal field structure (the poloidal flow structure in the hydrodynam- 
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ics) from the problem of particle acceleration and the structure of electric 
currents. The solution of the latter task for a given poloidal field can be 
obtained in terms of quite simple algebraic relations. It is important that 
such an approach can be straightforwardly generalized to flows in the vicin- 
ity of the rotating black holes, as the Kerr metric is also axially symmetric 
and stationary. 

On the other hand, it is much more difficult to find the two-dimensional 
poloidal magnetic field structure (the hydrodynamical flow structure) . First 
of all, this is due to the complex structure of the equation describing ax- 
isymmetric stationary flows. In the general case, it is a nonlinear equation 
of the mixed type, which changes from elliptical to hyperbolical at singu- 
lar surfaces and in addition contains integrals of motion in the form of free 
functions. Generally speaking, similar equations, which stem from the clas- 
sical Tricomi equation, have been discussed since the beginning of the last 
century in connection with transonic hydrodynamic Hows 3 . Later on, the 
equations describing axially symmetric stationary flows were called Grad- 
Shafranov equations after the authors who formulated in the late 1950s an 
equation of such a type in connection with controlled thermonuclear fu- 
sion 4 . This equation, however, was originally related to equilibrium static 
configurations only and required strong revision when it was generalized to 
the transonic case. The full version of such an equation was formulated by 
L.S. Soloviev in 1963 in the third volume of Problems of Plasma Theory 5 
and was well-known to physicists. However, as it often occurs, the full ver- 
sion of the Grad-Shafranov equation was little known in the astrophysical 
literature, so it was later 'rediscovered' scores of times 6 . 

As it turned out, the difficulty lay in the fact that the very formulation 
of the direct problem within the Grad-Shafranov approach proved to be 
nontrivial. For example, in the hydrodynamical limit, when there are only 
three integrals of motion, the problem requires four boundary conditions 
for the transonic flow regime. This implies that, for instance, two ther- 
modynamic functions and two velocity components should be specified at 
some surface. However, to determine the Bernoulli integral, which naturally 
should be known in order to solve the Grad-Shafranov equation, all three 
components of the velocity must be specified, which is impossible since the 
third velocity component itself is to be obtained from the solution. This is 
in fact one of the main difficulties of the approach under consideration. 

Nevertheless, several approaches exist that allow us to construct the 
analytical solution of direct problems within the framework of the Grad- 
Shafranov method. For example, this is possible when the exact solution 
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of this equation is known and we explore flows weakly diverging from the 
known one. Spherically symmetric accretion (ejection) of matter could 
be such an exact solution. As a result, the known structure of the flow 
in the zeroth approximation enables us to determine (with the required 
accuracy) both the location of singular surfaces and all the integrals of 
motion directly from boundary conditions, thus making it possible to solve 
the Grad-Shafranov equation within the direct formulation of a problem. 



1. Grad-Shafranov equation 

Let us consider the axisymmetric stationary plasma flow in the vicinity of 
a rotating black hole, i.e., in the Kerr metric 2 : 

ds 2 = -a 2 dt 2 + g ik (dx i + f3 l dt)(dx k + f3 k dt), (1) 

where 

a=£VA, /r = /9» = 0, p = - u = -™*, 

n 2 

"22 /r>\ 

9rr = 999 — P , g vv = W . (Z) 

Here a is the lapse function (gravitational red shift) vanishing on the hori- 
zon 



r s = M+ y/M 2 -a 2 , (3) 

w is the angular velocity of local nonrotating observers (the so-called Lense- 
Thirring angular velocity), and 

A = r 2 + a 2 - 2Mr, p 2 = r 2 + a 2 cos 2 9, 

E 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 9, w=-sm6. (4) 

p 

As usual, M and a are the black hole mass and angular momentum per unit 
mass (a = J/M) respectively. Here indices without hats denote components 
of vectors with respect to the coordinate basis d/dr, 3/ 89, and 8 /dip, and 
indices with hats correspond to their physical components. Finally, below 
we shall use the system of units with c = G = 1 . 

In what follows we shall also use the 3 + 1 split language 2 . Within 
this approach, the physical quantities are expressed in terms of three- 
dimensional vectors which would be measured by observers moving around 
the rotating black hole with the angular velocity ui (so-called ZAMOs - zero 
angular momentum observers). The convenience of the 3 + 1 split language 
is connected with the fact that the representation of many expressions has 
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B = - (5) 



the same form as in the flat space. On the other hand, all thermodynamic 
quantities are determined in the comoving reference frame. 

Now, we shall demonstrate how the five 'integrals of motions', which 
are constant at the magnetic surfaces, can be derived in the general case 
of axisymmetric stationary flows. It is convenient to introduce the scalar 
function \I/(r, 9) which has a meaning of magnetic flux. As a consequence, 
the magnetic field is defined in the following way: 

VjE x _ 2J_ 
2ttw azu 

where I(r, 9) is the total electric current inside the region ^ < ^(r, 9). 

As usual, we assume that the magnetosphere contains sufficient amount 
of plasma to satisfy the freezing-in condition which, using the 3 + 1 split 
language, preserves the form E + v x B = 0. On the other hand, the 
stationarity (as well as the condition for zero longitudinal electric field) 
implies that the field E can be written as 

E = -^P^V*. (6) 
27ra v ; 

By substituting relation (6) into the Maxwell equations, it is easy to verify 

that the condition B • VSlp = is satisfied, i.e. that f2p must be constant 

at the magnetic surfaces (Ferraro's isorotation law): 

fi F = fi F (*). (7) 

Next, the Maxwell equation V • B = 0, the continuity equation, and the 
freezing-in condition allow us to write the four-velocity of matter u in the 
form 

u= -!LB + "f(n F -Lu)— e # , (8) 
an a 

where 7 = — v 2 is the Lorentz factor of matter (measured by 

ZAMOs), and the quantity 77 is the particle flux to magnetic flux ratio. 

Due to the relationship V • (nH p ) = 0, n must be constant at the magnetic 

surfaces ^(r, 9) — const as well, i.e., 

77 = ??(*). (9) 

The next two integrals of motions result from our assumption that the 
flow is axisymmetric and stationary. This yields the conservation law of 
the energy E and the z-component of angular momentum L z : 

E = E{^) = + fin(aj + uju^); (10) 
L = L(V) = +nnzuu l f,, (11) 

Z7T 
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where /i = (p m + P) jn is the relativistic enthalpy (p m is the internal energy 
density, and P is the pressure). Finally, in the axially symmetric case the 
isentropy condition yields 

* = «(*), (12) 

so that the entropy per particle, s(v&), is the fifth integral of motion. 

The five integrals of motions Qf(*)j »?(*)) £(*), and £(*), as 

well as the poloidal magnetic field B p , allow us to find the toroidal magnetic 
field B(p and all other plasma parameters: 

/ _a 2 L-(n F -Lu)m 2 (E-ojL)^ 



2ir a 2 - (0 F - lu) 2 w 2 - M 2 



(13) 



1 a 2 (E-n F L)~M 2 (E-ujL) ^ 
7 ~ a 2 - (0 F - u) 2 w 2 - M 2 ' ( ' 



1 (E - n F L)(n F - uj)zu 2 - LM 2 
voyLr\ a 2 - (0 F - lu) 2 vj 2 - M 2 



(15) 



where 



M 2 = (16) 



It is easy to see that M 2 is proportional (with the factor of a 2 ) to the 
Mach number squared of the poloidal velocity u p with respect to the Alfven 
velocity ua — B p / ^/Aimp, i.e., M 2 = a 2 u 2 /u\. 

Since n = fi(n, s), definition (16) allows us to express the concentration 
n (and hence the specific enthalpy fi) as a function of r], s, and M 2 . This 
means that along with the five integrals of motion, the expressions for I, 7, 
and U(p depend only on one additional quantity, namely the Mach number 
M. To determine the Mach number M, it is necessary to use the obvious 
relation -f 2 — u 2 = 1, which, owing to equations (14) and (15), can be 
rewritten in the form 

K 1 M 4 (W) 2 2 2 2 

^ = 64^— b^ +aW ' (17) 

where 

A = a 2 - (0 F - cufrn 2 - M 2 (18) 



and 



K = a 2 w 2 (E - F L) 2 [a 2 - (0 F - lo) 2 w 2 - 2M 2 ] 

+M i [w 2 {E-uL) 2 -a 2 L 2 } . (19) 
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As for the Grad-Shafranov equation itself, i.e., the equilibrium equation 
for magnetic field lines, it can be written in the form 



d* 
(20) 



647T 4 1 d (G\ ~ ld V 1R , ^ 



+ o n TTTTo TTT" ~T ~ 167T>n- — f - 167T nT— — = 0, 

a 2 m 2 2M 2 d^\AJ p 77 d* d* 



where 

G = a 2 n7 2 (£ - F L) 2 + a 2 M 2 L 2 - M 2 zu 2 (E - loL) 2 , (21) 

and the derivative d/d^> acts on the integrals of motion only. Finally, 
expressing in Eq. (20) the terms V&M 2 according to Eq. (19) we obtain 7 

Iv ( 1 v k A 1 1 v^-Vtt-ViVfctt 

a H«^ 2 / a 2 w 2 (Vi<) 2 £> 

+ V ^ V ^ £ -^V' fe FV^ + ^^(W) 2 ^ (22) 

a 2 w 2 a 2 w 2 (V*) 2 2D a 2 v ' d* v ' 

64tt 4 1 9 /G\ ,„ , 1 dry „ , Hs „ 

+^^7^7?Trr t - 16tt 3 ^-tt - 167r 3 nT— = 0. 
a 2 w 2 2M 2 9* V A / r\&$> d* 

Here 

^ = A + ^14_^A^L (23) 

M 2 M 2 B 2 u 2 M 2 1 - c 2 ' k J 

P P a 

^ 64tt 4 X 64tt 4 2 2 2 2 
= Ivi 4 " A 2 M^" (24) 
and the gradient Vj-. denotes the action of under the condition that M 
is fixed. Let us stress that in equation (22) the pressure P, the temperature 
T, sound velocity c s , and the specific enthalpy /x are to be expressed via an 
equation of state in terms of the entropy s(\P) and the square of the Mach 
number M 2 . In turn, the quantity M 2 is to be considered as a function of 
(V*) 2 and the integrals of motion, 

M 2 = M 2 [(V*) 2 ,S(*),L(*),77(*),r! F (*), S (*)] . (25) 

The latter relation is the implicit form of Eq. (17). The stream equation (22) 
coupled with definitions (5) - (11) is the desired equation for the poloidal 
field which contains only the magnetic flux 'J and the five integrals of 
motion fi F (*), »?(*), and L(\I>) depending on it. 
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Equation (22) is a second-order equation linear with respect to the high- 
est derivatives. It changes its type from elliptical to hyperbolical at singular 
surfaces where the poloidal velocity of matter becomes equal to either fast 
or slow magnetosonic velocity (D = 0), or to the cusp velocity (D = — 1). 
Although at the Alfvenic surface, A = 0, the type of equation does not 
change, the Alfvenic surface does represent a singular surface of the Grad- 
Shafranov equation because a regularity condition must be satisfied there. 



2. Examples 

Bondi-Hoyle accretion. As a first example, we consider the hydrody- 
namic accretion onto a moving back hole (the Bondi-Hoyle accretion) , which 
is one of the classical problems of modern astrophysics 2 . First of all, let 
us formulate the hydrodynamical limit of the Grad-Shafranov equation, 
where we can neglect the electromagnetic field contribution. In this case, 
it is convenient to introduce a new potential $(\&) satisfying the condition 
rj(^) = d^/d*. Using definition (8) we obtain 

cmu p = — !— (V$ x e«A (26) 
Zttzu 

Surfaces <I>(r, 9) — const define the streamlines of matter. 

In the hydrodynamic limit, there are only three integrals of motion. 
These are the energy flux and the z-componcnt of the angular momentum: 

E(§) = (i(a~f + zuluu ); (27) 
L($) = [ivou^,, (28) 

as well as the entropy s = s($). Now the algebraic Bernoulli equation (17) 
takes the form 



a 2 ^ M 4 



where the 'Mach number' squared M 2 is defined as M 2 = 4irn/n. Then 
the Grad-Shafranov equation (20) can be rewritten in the form 7 

/ d$ 




zu (E - uiL) — - u— - a z L— 



(30) 



where now 



° = -i + ^- < 31 > 
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As we see, equation (30) contains only one singular surface, i.e. the sonic 
surface, determined from the condition D — 0. 

To construct the solution corresponding to the Bondi-Hoyle accretion, 
it is possible to seek the solution of the Grad-Shafranov equation for the 
flux function <I>(r, 9) in the form of a small perturbation of the spherically 
symmetric flow in the reference frame moving with the black hole 

$(r,6>) = $ o [l-cos0 + £i/(r,0)]. (32) 

Here we introduce a small parameter 

e x = ^ (33) 

which defines the ratio of the black hole velocity to the velocity of sound 
at infinity. For a nonmoving gravity center we return to the spherically 
symmetric flow. 

As Grad-Shafranov equation (30) contains three invariants, it is neces- 
sary to specify four boundary conditions, say 

(1) w P ,oo = const, 

(2) v v = (and hence L = 0), 

(3) Sqo = const, 

(4) B 00 =4,/(T-1). 

In the last relation we neglect the terms ~ e\. 

As a result, the Grad-Shafranov equation can be linearized: 



-e^D^ - -L { D + 1) sin.- ^—^ j + s^N^ = 0, (34) 



where 

2 [i 1 M 

It is extremely important that according to (26) and (31), 

D + l= " V 2 2 -t^, (36) 

so the factor a 2 enters every term of equation (34). Hence, equation (34) has 
no singularity at the horizon. In particular, it means that it is not necessary 
to specify any boundary conditions for r — r s . It is not surprising because 
the horizon corresponds to the supersonic region which cannot affect the 
subsonic flow. 
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On the other hand, we see that all the terms contain the small value 
£\. Hence, the functions D, c s , etc. can be taken from the zeroth solution. 
As for the spherically symmetric flow, the functions D, c s , etc. do not 
depend on 9, and the solution of equation (34) can be expanded in eigen 
functions of the operator sin 9 d/ 88(1/ sin9 ■ d/89). Thus, the solution can 
be presented in the form 

oo 

f(r,e)=J29m(r)Q m (0), (37) 

m=0 

the equations for the radial functions g m (r) being 

2 7-1°- 9m . 2 i\t °-5m i , \ A* c s r\ too\ 

r D—rr+ r N r—, \-m(m+l)— — g m = 0. 38 

Here Qq = 1 — cos 9, Q\ = sin 2 #, Q 2 = sin 2 9 cos 9, ... are the eigen 
functions of the angular operator. 

As to the boundary conditions, they can be formulated as follows: 

(1) No singularity on the sonic surface (where N r = 0, D = 0), 

9m (r,) = 0. (39) 

(2) The homogeneous flow $ = TmocVoor 2 sin 2 9 at infinity which gives 

9i^7^ 2' 52,53,- -- = 0. (40) 

As a result, the complete solution can be presented in the form 

$(r,0) = $ o [l-cos0 + ei 5 i(r)sin 2 6»)], (41) 

where the radial function g\{r) is the solution of the ordinary differential 
equation (38) for m = 1 with the boundary conditions (39) and (40). 

This means that we have constructed the analytical solution of the prob- 
lem, i.e., obtained the full description of the flow structure. For example, 
the sonic surface has now the nonsphcrical form 



r*(0) = r* 



T + 1 

1 + £i ( 5 _ 3r ) fecostf 



(42) 



where the numerical coefficient k 2 = r*^(r*) is expressed through the 
derivative of the radial function gi(r) at the sonic point. As shown in 
Fig. f , the analytical solution fully agrees with the numerical calculations 10 
in spite of the parameter E\ = 0.6 here being quite large. Here r = 4/3, 
and the numbers alongside the curves denote values $/$o; the dashed lines 
show the streamlines and the form of the sonic surface obtained numerically. 
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Figure 1. Bondi-Hoyle accretion onto moving black hole. 

As can be easily seen, outside the capture radius 

R c sa e~ 1/2 r* (43) 

our main assumption, i.e., the smallness of the deviation from the spheri- 
cally symmetric flow, is not valid. Nevertheless, the solution found remains 
correct. This remarkable property is due to the Grad-Shafranov equation 
becoming linear for constant concentration n. But as we learn from the 
spherically symmetric Bondi accretion, at large distances r 3> from the 
sonic surface the density of the accreting matter is virtually constant. Ac- 
cordingly, the concentration is constant for a homogeneous flow as well. As 
a result, under the condition that R c ^> r*, which holds true for e% <C 1, 
near and beyond the capture radius (where the perturbation ~ £i<7i(r) be- 
comes comparable to unity) equation (30) becomes linear. So that the sum 
of the two solutions, homogeneous and spherically symmetric ones, is also 
a solution. 

Thin transonic disk. As a next example, we consider the internal 
two-dimensional structure of a thin accretion disk. Here, for simplicity we 
consider the case of a nonrotating (Schwarzschild) black hole 11 . We re- 
call that according to the standard model 12 the accreting matter forms an 
equilibrium disk rotating around the gravitational center with the Keple- 

1/2 

rian velocity Vfc(r) — (GM/r) . The disk will be thin provided that its 
temperature is sufficiently small (c s <C Vjt) since the vertical balance of the 
gravity force and the pressure gradient implies that 

H^r—. (44) 

Vk 
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The General Relativity effects result in two important properties: the ab- 
sence of stable circular orbits for r < ro = 3r g and the transonic regime 
of accretion. The first point means that the accreting matter passing the 
marginally stable orbit approaches the black hole horizon sufficiently fast, 
namely, in the dynamical time Td ~ [u r (ro)/c] _1 / 3 r g /c. ft is important 
that such a flow is realized in the absence of viscosity. The second state- 
ment results from the fact that up to the marginally stable orbit the flow 
is subsonic while at the horizon the flow is to be supersonic. 

It is necessary to stress that the existence of the small parameter e 2 = 
uo/co <C 1, where c is the sound velocity and uq is the gas radial velocity on 
the marginally stable orbit, comes from the relation w r /«K ~ c^ssCg/v^ for 
the radial velocity in the accretion disk 1 . In the vicinity of the marginally 
stable orbit this estimate is apparently inapplicable. Nevertheless, below 
we shall consider the parameter Si to be small because the presence of 
a small parameter allows us to investigate the flow structure analytically. 
In addition, the small parameter makes the effect under discussion more 
visible. 

Up to now in the majority of works devoted to thin accretion disks the 
procedure of vertical averaging was used, where the vertical four-velocity 
Ug was assumed to be zero 13 . As a result, the vertical component of the 
dynamic force nu h Vb(^u a ) in the Euler equation was postulated to be inim- 
portant up to horizon. For this reason the disk thickness was determined by 
the pressure gradient even in the supersonic region near the black hole 14 . 
Here I am going to demonstrate that the assumption ug = is not correct. 
As for the Bondi accretion, the dynamic force is to be important in the 
vicinity of the sonic surface. 

Figure 2 shows the structure of a thin accretion disk after passing the 
marginally stable orbit r = 3r g obtained by solving equation (30) numeri- 
cally for Co = 10~ 2 , Mo = 10~ 5 , i.e., in the presence of a small parameter 
£2 = Mo/co <C 1. The solid lines correspond to the range of parameters 
u p/ c o < 0-2- The dashed lines indicate the extrapolation of the solution to 
the sonic-surface region. In the vicinity of the sonic surface the flow has 
the form of the standard nozzle. 

As we see, the flow structure near the sonic surface is far from being 
radial. The appearance of the narrow waist has a simple physical meaning. 
Indeed, the density remains almost constant for subsonic flow while the 
radial velocity increases from u to c* <~ c , i.e., for £2^1 changes over 
several orders of magnitude. As a result, the disk thickness H should change 
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Figure 2. Nozzle shape of a thin transonic disk near marginally stable orbit. 

in the same proportion owing to the continuity equation as well, 

ff(r.) « -H(3r g ). (45) 
Co 

Here it is extremely important that both components of the dynamical force 
become comparable with the pressure gradient near the sonic surface: 

~ Uf ~ ~ o ' 

r at) or /i Mq r 

where the angle 6 is counted off from the equatorial plane. 

In other words, if there appears a nonzero vertical velocity component, 
the dynamical term (vV)v cannot be neglected in the vertical force bal- 
ance near the sonic surface 11 . It is clear that this property remains valid 
for arbitrary radial velocity of the flow, i.e., even when the transverse con- 
traction of the disk is not so pronounced. Taking dynamical forces into 
account causes two additional degrees of freedom to appear, which relate 
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to the higher derivatives in the Grad-Shafranov equation. This also leads 
to extra conclusions independent of the value of £2- In thin accretion disks, 
the critical condition at the sonic surface does not fix the accretion rate any 
more; it determines the bending of the streamlines near the sonic surface. 

Finally, the inclusion of the vertical velocity inevitably leads to the 
appearance of a small longitudinal scale Sru w if* in the vicinity of the 
sonic surface, which for a thin disk proves to be much smaller than the 
distance to the black hole for any value of the parameter uo/cq. In the 
standard one-dimensional approach, this scale does not emerge. As for the 
supersonic region (and, in particular, the region near the horizon), the disk 
thickness here will be determined not by the pressure gradient, but by the 
form of ballistic trajectories. 

The Blandford-Znajek process. In conclusion, we discuss the energy 
loss of a rotating black hole embedded in an external magnetic field - the 
so-called Blandford-Znajek process 15 . This process is considered to be 
the most preferential mechanism of energy release in active galactic nuclei, 
microquasars, and even cosmological gamma-ray bursters 16 . Its main idea 
is based on the analogy with the energy transfer in the internal regions of 
radio pulsar magnetospheres. Indeed, let us suppose that there is a regular 
external magnetic field in the vicinity of a rotating black hole, and that the 
electric current I flows along magnetic field lines. Then, the electric field 
E, which is induced by plasma rotating with the angular velocity and 
the toroidal magnetic field B v , which is due to the longitudinal current I, 
generate the electromagnetic energy flux (the Poynting vector flux) carrying 
the energy away along the magnetic field lines. 

Of course, by definition, general relativity effects are important near 
the black hole. Consequently, it is not obvious that the pulsar analogy 
can be useful in all cases. For example, in radio pulsars, the braking of 
neutron stars results from the Ampere force acting on the star surface. This 
force results from the surface currents shorting the electric currents flowing 
in the pulsar magnetosphere 17 . In the case of black holes such currents 
cannot lead to deceleration as the event horizon is not a physically preferred 
surface, though surface currents themselves can be formally introduced in 
the framework of the so-called membrane paradigm 2 . 

Indeed, let us consider the well-known condition at the horizon (the 
absence of infinite electromagnetic fields in the reference frame comoving 
with freely falling observer) E' g — ► (Eg + B v )/a < 00. It means that 
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Eg + B v — > which can be rewritten in the form 

r 2 + a 2 /Hv|/\ 
47TJW = [0 H - F (^)] % sing — . (47) 

r; + cos^ 6* V at* / „ 

Here a = J/M is the rotation parameter. This relation was used 15 as the 
'boundary condition on the horizon'. But it is clear that the horizon is not 
causally connected with the outer space. For this reason the conclusion was 
made that there is no energy flux along magnetic field lines passing through 
the black hole horizon 18 , and hence a black hole cannot work as a unipolar 
inductor extracting the rotation energy by electromagnetic stresses. 

However, a recent more accurate analysis 19 (in which, in fact, the first 
solution of the Grad-Shafranov equation for nonzero particle mass in the 
Kerr metric was obtained) indicates that the braking torque acts in the 
plasma generation region above the black hole horizon. Such a torque 
appears due to the action of long-range gravitomagnetic forces which pene- 
trate into the regions causally connected with the outer space. As a result, 
it was shown that for a finite mass of particles in the very vicinity of the 
horizon, there is a hyperbolic region of the Grad-Shafranov equation which 
is altogether absent in the force-free approximation. Hence, the full ver- 
sion of the Grad-Shafranov equation needs no boundary condition on the 
horizon. Actually, this property was already demonstrated above by the 
example of the Bondi-Hoyle accretion. Thus, in this case it is impossible 
to consider relation (47) as a boundary condition. Relation (47) is auto- 
matically true for any solution of the Grad-Shafranov equation which can 
be extended up to the horizon. 

On the other hand, in the force-free approximation, when the Grad- 
Shafranov equation remains elliptical up to the black hole horizon, extra 
condition (47) is to be included into consideration 20 . But as was demon- 
strated 9 , this condition is actually the manifestation of the critical con- 
dition on the fast magnetosonic surface. Hence, in reality this condition 
is specified on the surface which does not coincide with the event horizon, 
and, hence, is in the causal connection with the outer magnetosphere. 

Thus, the Blandford-Znajck mechanism of electromagnetic energy ex- 
traction from rotating black holes faces no causality problem. As in the 
pulsar magnetosphere, if there is enough secondary plasma to screen the 
longitudinal electric field, its charge density and electric currents produce 
the flux of electromagnetic energy propagating from the central star to in- 
finity. For the same reason, a rotating black hole embedded into an external 
magnetic field works as a unipolar inductor extracting its energy of rotation 
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by the flux of the electromagnetic energy. As a result, the energy loss can 
be evaluated as Wtot ~ Wbz, where 



W BZ 



fi F (Sl H - F ) / a \ 2 2 2 



ft 2 VmV 



It is easy to check that for the extreme rotation of a black hole (a = 1) and 
for B = S Ed d « 10 4 (M/10 9 M Q )- 1 /2 Q, the energy loss W B z (48) coincides 
with the Eddington luminosity. 




a) b) 

Figure 3. Magnetospheres of rotating black hole (a) and radio pulsar (b). 



It should be emphasized, however, that as follows from equation (48), 
the rate of the energy release needed to explain the luminosity of ac- 
tive galactic nuclei can be achieved only for the extreme black hole mass 
~ 10 9 M Q , the extreme magnetic field B ~ i?Edd, and the extreme angu- 
lar velocity a ~ M. Therefore, some papers have recently appeared in 
which the efficiency of the Blandford-Znajek process in real astrophysical 
conditions was questioned 21 . In particular, it was pointed out that for 
rapid rotation the Wald solution for the vacuum magnetosphere leads to 
the magnetic field being pushed out from the horizon into the ergosphere 2 , 
which causes the appearance of the additional factor (1 — a 2 /M 2 ) — > in 
expression (48). 

But as shown in Fig. 3a, in the black hole magnetosphere filled with 
plasma, all magnetic field lines crossing the surface of the internal light 
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cylinder a 2 = (flp — lo) 2 w 2 + M 2 are to cross the black hole horizon as 
well, that is why to an order of magnitude the energy loss for extremely 
rotating black holes coincides with the loss given by equation (48). Indeed, 
here the situation is to be fully analogous to the pulsar magnetosphere 
where the field lines passing the external light cylinder do not intersect the 
equatorial plane (Fig. 3b). The magnetic field structure shown in Fig. 3a 
was recently obtained numerically 22 . 

3. Conclusion 

In some simple cases the Grad-Shafranov equation allows us to construct 
the exact solution to the problem. In particular, this approach is very 
useful in studying the analytical properties of transonic flows and in deter- 
mining the required number of boundary conditions. On the other hand, 
in the general case no consistent procedure exists regarding the construc- 
tion of the solution within the Grad-Shafranov approach. The point is 
that the location of singular surfaces, at which critical conditions should 
be formulated, is not known beforehand and itself must be found from the 
solution to the problem. Moreover, it is impossible to generalize this ap- 
proach to the case of nonideal, non-axially symmetric and nonsteady flows. 
So it is not surprising that most investigators, who are in the first place 
interested in astrophysical applications, have focused on a totally different 
class of equations, namely on time relaxation problems, that can only be 
solved numerically 23 . Nevertheless, it is clear that the key physical results 
obtained using the Grad-Shafranov approach are independent of the com- 
puting method. For this reason one can hope that the results presented 
above can be useful for everyone working in this field. 
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